Delocalisation transition of a rough adsorption-reaction interface 
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We introduce a new kinetic interface model suitable for simulating adsorption-reaction processes 
which take place preferentially at surface defects such as steps and vacancies. As the average interface 
velocity is taken to zero, the self- afBne interface with Kardar-Parisi- Zhang like scaling behaviour 
undergoes a derealization transition with critical exponents that fall into a novel universality class. 
As the critical point is approached, the interface becomes a multi-valued, multiply connected self- 
similar fractal set. The scaling behaviour and critical exponents of the relevant correlation functions 
are determined from Monte Carlo simulations and scaling arguments. 
PACS Numbers: 68.35.Ct, 82.65. Jv, 05.70.Ln, 68.35.Rh 



Kinetically roughened interfaces display a rich phe- 
nomenology, have deep connections with fields as diverse 
as self-organized criticality, spin-glasses and complex pat- 
tern formation, and lend themselves to modelling various 
systems with practical applications, ranging from het- 
erogenous catalysis to geomorphology The huge 
amount of numerical and analytical effort that has re- 
cently been invested in them has revealed that they obey 
universal scaling relations, which fall into one of a few 
universality classes. In this letter we would like to present 
a kinetic interface model which exhibits an anisotropic to 
isotropic phase transition with novel scaling behaviour at 
the delocalisation critical point. 

Reaction fronts formed by A + B — > reactions in het- 
erogeneous systems where the reaction takes place on a 
two dimensional substrate || are often confined to a nar- 
row "reactive zone" especially if the reactants are either 
initially segregated or become segregated due to reaction 
kinetics [ 5^ 9j. Our present model is motivated by recent 
findings 1 1C |ll[] of high reaction rates and strong bonding 
at surface defects like steps and vacancies in studies of 
heterogeneous catalysis, a burgeoning new field in surface 
science. 

We consider an idealised surface with only one step, 
terminating a terrace made up of A particles (Figjl). 
The surface is exposed to two kinds of incoming parti- 
cles, A and B, which are allowed to adsorb at first con- 
tact, and only on sites adjacent to the step, which we 
will call "interface sites." The adsorption of A particles 
makes the interface advance. The adsorbing B particles, 
on the other hand, immediately react with an A neigh- 
bor to form a product which leaves the surface. This eats 
into the step, making the interface recede. We investi- 
gate the effect of changing the rate of injection [Q of 
the two reactants. We do not allow any reactions to take 
place with the substrate atoms. We assume, for simplic- 
ity, that the temperature is low enough so that no surface 
restructuring occurs; the bonding to the interface sites is 
sufficiently strong [Of for diffusion along the interface to 



be prohibited. The kinetics is therefore driven by the ad- 
sorption and reaction steps and not by the transport of 
the reactants. 




FIG. 1. A terrace of A particles, with the interface sites, 
indicated by dots, neighboring the step. 

The model is defined on an infinite strip of width L, 
on which we impose periodic boundary conditions. The 
interface is initially a perfectly straight line located at 
h = 0. The system is driven weakly so that at any instant 
only one particle of either A or B type, with probabilities 
Pa or pb = 1 — Pa, impinges on the interface. As the in- 
terface moves with a mean velocity equal to e = pa — 0.5, 
it roughens, and becomes multiply connected, shedding 
"islands" or "lakes" in its wake. As |e| — ► 0, the growth 
direction is completely delocalised, the width of the in- 
terfacial region keeps on growing indefinitely, and the in- 
terface breaks up into an isotropic fractal (see Fig.|J). For 
finite L, there may exist more than one spanning string 
of interface sites at e = ; this phenomenon is similar to 
the formation of Liesegang bands ]l4|| ). It is the purpose 
of this letter to understand the nature of this derealiza- 
tion transition, to describe the crossover behaviour and 
to characterise the self-similar reactive region formed as 
|e|-0. 

For many interface problems with a well defined growth 
direction, such as the Eden |l5[ model, or the Edwards- 
Wilkinson model fllql, where the interface can be de- 
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scribed with a single- valued, self-affine curve. The scal- 
ing behaviour of the interface width may be conveniently 
summarized by the scaling form 17 



l/z\ 



(1) 



where g(u) ~ const for u < 1 and ~ u x for u 3> l;z is 
the dynamical critical exponent, and (3 = xl z ■ Kardar, 
Parisi and Zhang |L8) have found the values z = 3/2, 
X = 1/2 and f3 = 1/3 for the stochastic differential equa- 
tions describing Eden growth in d = 1 + 1. This set of crit- 
ical exponents characterizes a wide range of anisotropic 
growth phenomena with annealed noise [jjj, and where 
the local velocity of the interface increases with the slope. 
In the limit that the velocity goes to zero or is indepen- 
dent of the slope ||, one gets the Edwards- Wilkinson 
model |l(| , which is exactly solvable in d = 1 + 1 dimen- 
sions and falls into another universality class, character- 
ized by z = 2, x = 1/2 and (i = 1/4. 




FIG. 2. The interface at e = 0, for early and later stages 
of growth. The configurations for early times resemble the 
surface at larger values of e. 

Since our interface is typically multivalued we define 
the width function within an interval of size £ as, 
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where hi is the height of the ith interfacial site, i — 
1, . . . , N(£), and h(£) the mean position of the interface. 
We find j3 = 1/3 for early times. In the limit of pa = 
or pb — 0, i.e., for |e| — > 0.5, our model is equivalent 
to Eden |ll| growth, and indeed, along the singly con- 
nected part of the interface x — 1/2 in the steady state. 
However, effective roughness exponent (Xeff) goes contin- 
uously to zero as |e| — > as shown in FigJ|. 
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FIG. 3. The effective roughness exponent Xcff , for discon- 
nected parts of the interface included in (o) and excluded from 
(•) the analysis; L < 1024. The error bars are comparable to 
the size of the symbols. 

The reason for this is that as we decrease |e|, the sur- 
face becomes highly convoluted, with islands (or lakes) of 
all sizes, and therefore increasingly multivalued. A com- 
peting length scale emerges in the system, the "thick- 
ness" of the interface, which we can measure by the vari- 
ance of the height, 
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with rii being the number of interfacial sites {hij}, above 
any point i along the horizontal axis. The thickness 
obeys [^9| a skewed-Gaussian distribution. The average 
yi, =< y >l and the second moment of this distribu- 
tion both diverge as |e| — > with a critical exponent 
v = 0.55 ± 0.05 ~ 1/2, as - . The scaling form for 
Vl is, 
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where G(v) ~ v for v < 1 while for v > 1, G(v) ~ const, 
with = 1/2 and the (longitudinal) dynamical critical 
exponent £ obeys £ = 1//3. The thickness of just the 
singly connected part does not diverge as |e| — > 0, so 
that the disconnected parts make up almost all of the 
interfacial region at the derealization transition jl9| . 

Normalizing w(£) in Eq.(||) by y 1 ^ 1 yields a collapse of 
the data for all e as can be seen from Fig.||. We believe 
that the small deviation of the power of t/l from unity 
is due to insufficient statistics for y^ which converges ex- 
tremely slowly as |e| — > and may safely be neglected, 
and we conclude 
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const. £ 1/2 < y L . 
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Thus, Xeff goes to zero as the self-affine excursions of the 
interface are blurred by the thickness of the interface as 
yL becomes greater than i 1 ' 2 . 
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FIG. 4. The width normalized by a power of the thickness 
t/L very near unity, displays KPZ behaviour for l 1 ^ 2 /yL > 1. 
The different curves correspond to data taken at e = 0.001 , 
0.003, 0.005, 0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5. The deviations 
from the smooth collapse are due to I becoming comparable 
to the system size L = 1024. 

If one considers coarse-grained width functions, either 
by taking the average height at any given point or the 
maximum pC| ] height for e > (minimum for e < 0), 
one finds that they obey the scaling form (|]) with KPZ 
exponents for £ <C L. 

To get the local scaling picture, we focus on a single 
spanning string in the interface and consider 



C x (l) = 
C h {l) = 
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where both r and I are the ( "chemical" ) length measured 
along the string, and x and h are Cartesian coordinates 
of the interface site. The scaling relations we have found 
from Fig|5| for these quantities are given below. In the 
transient regime (I ^> t 1 / 2 ), 



I 

I ft* 



VL « I 112 
VL » I 112 



and 



(7) 



(8) 



where if) = 1/6. 

Note the horizontal projection of a segment of fixed 
"chemical length" decreases with t as the surface crum- 
ples with time in the critical (|e| — > 0) region. In the 
steady state, (I <C t 1 / 2 ), 
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We see that for 3> I 1 / 2 , the interfacial region becomes 
isotropic, with C x ~ Ch- This regime is characterized by 
an "isotropic" roughness exponent Xisot = 2/3. In the 
opposite limit, Ch ~ Cx' 2 , as expected for the self-affine 
Eden surface. This crossover is clearly seen in FigJ|. 
In accordance with the above observations we propose 
the following scaling functions for the whole range of e. 
Defining u = l/t 1 ^ and s = i/l/I x , we have, 



C h ~a h (e) t & f(u,s) 
C x ~ a x (e) t 1/z g(u,s) 



where 
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u > 1 

u < 1, s < 1 
u < 1, s > 1 



s < 1 
u > 1, s > 1 
it < 1, s > 1. 
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where x, z, and (3 have their KPZ values. The ampli- 
tudes are defined as a x (e) = (a + lei 1 / 6 ) and ah — 
l/(a _1 + lei 1 / 6 ), where a is some constant. 
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FIG. 5. a) C h and (7* for e = 0.5 for I = 2, L/2, and 
different times, b) Ch in the isotropic region, e = 10 -4 . 

From these scaling forms and Eq.(||) we see that the 
new critical exponents obey the relationships, 
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FIG. 6. The vertical projection Ch of sections of the singly 
connected part, plotted as a function of the horizontal pro- 
jection C x . 

From Eqs.fl|lO|) we see that the graph dimension of 
a singly connected part in the isotropic regime is D g = 
\/ X lsot 

roughness exponent via D g 



Since in two dimensions, D g is related to the 



2 — x, for x — 1/2 we get 



X isot = 2/3. The scaling relation <\l§ yields 13= 1/2, 
from and (|l|) it follows that ip = 1/6 and C = 2. 
The fractal dimension of the self-similar set of interface 
sites JilJ within a band of width ~ \e\~ u is found, 
from boxcounting to be Dj = 1.85±0.05 for length scales 

t<vh 

In conclusion, we have presented an absorption- 
reaction model where the interface undergoes a delocal- 
isation transition at the point where the mean velocity 
of the interface goes to zero. Although it might be con- 
jectured Q that as the velocity of the interface vanishes, 
the scaling behaviour should cross over to the Edwards- 
Wilkinson universality class, this is not the case here. It 
has previously been observed [^,|lj, that the presence 
of overhangs, islands and inclusions may cause the small- 
scale structure of the interface to crossover form being 
self-affine to self-similar while the large scale behaviour 
remains self-affine. In the present model, this crossover 
is driven by a competing length scale, the thickness of 
the interface, which diverges at the critical point as \e\~ v 
with v = 1/2. In the critical region, the interface is char- 
acterized by new set of exponents x lsot = 2/3, $ = 1/2, 
£ = 2, rp = 1/6, and the fractal dimension Di = 1.85. 
Except for v and D/, these exponents may be obtained 
from the KPZ exponents via scaling relations. 

It should finally be mentioned that the reaction re- 
gion can be described by stochastic differential equations 
of the multiplicative noise type with a single component 
field, since in our model interface sites cannot be created 
spontaneously either in the bulk or the vacant region. 



A field-theoretic renormalization group computation a la 
Tu, Grinstein and Munoz p2fl is presently under way to 
obtain the values of the critical exponents. 
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